C:\Docuinents and SettingsX jhood\Local ... Files\OLK8\ReadOnly-0-Rev-13-ColorLearn . c 



1 



//////////////////////////////////////////////////////////////////////////////// 
II 

II Title : ColorLearn.c 

// Author : Siming Lin 
// Date : flBHi 

// Copyright : National Instruments iHV* All Rights Reserved. 
// Access : Company Confidential 

// Purpose : Implements the ColorHistogram function. 
// 

//////////////////////////////////////////////////////////////////////////////// 



// Include Directives 
# include PRECOMP 

#include "Analysis/Color/CommonColor .h" 
#include <stdlib.h> 
# include <math . h> 

//////////////////////////////////////////////////////////////////////////////// 



// 

// ColorLearn 
// 

// Description: 

// Extract the color feature from the input color image. 

// Feature: 

// the estimated color spectrum 

// Parameters: 

/ / pSrcRef - source image 

// pCoordPtr - optional rectangle 

// pNbColor - level of color representaion : 0, 1, and 2. 

// pSatThreshold - Threshold for distinguishing two colors with the same hue value 

// 

// pColSpectrumRef - calculated color spectrum feature 

// 

// Return Value: 
// Error code 



//////////////////////////////////////////////////////////////////////////////// 

// algorithms for color feature extraction 

enum ColorFeatureSet { 
kColSpectrum = 0, 
kColSignature, 
kColHis Intersection 
}; 



enum ColorComplexity { 
kComplexityLow = 0, 
kComplexityMedium, 
kComplexityHigh 
} ; 

#define BWHueThreshold 10 // this one may need be changed to adapt to light and hue 
#define BlkLgtThreshold 10 // need more experiment with IMAQ color board here ! 
#define BWLgtMiddle 12 8 
#define NbColorBinBase 7 

#define OffsetColor 18 // center of the red color for LOW complexty 



GRLIBError ColorLearn (const GRImage* pSrcIMRef, const ROI* pRoi, int pNbColor, int 
pSatThreshold, ArraylD* pColSpectrumRef) { 
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GRImage* IMasklmagePtr = NULL; 
MaskPart IMaskPart; 



const Pix8* iMaskPixPtr; 
const Pix8* IMaskLinePtr; 



const GRImage* 



ISrcIMPtr; 
iSrcPixPtr; 



const PixRGB 



*lSrcPixRGBPtr, *lBuf SrcRGBPtr; 



const PixHSL 



*lSrcPixHSLPtr, *lBuf SrcHSLPtr; 



double 



*lColSpectruniAlDWPtr; 



int 



i, ISizeXLW, ISizeYLW, ISrcLineWidth, IColLW, ILineLW 



double 



llmageAreaDW, IHueStepDW; 



int 



INbColFeature, INbHueBin, IHuelndex, lNumColorEntry=0 



long IRHHVal, IGSSVal, IBLVVal; 

long Ihue, Isat, llgt, IBlkThreshold; 

Byte Ireplace = 191, offset=0; 

long Icoring = 0; 

double h,s,l; 

long IMskLineWidth; 

int lpixelsOutsideMask=0; 

GRLIBError error = ERR SUCCESS; 



// lookup table for function f (x) =exp {-0 . 025*x) , x=0 - 200 



static double expLookup[] = { 





1. 


.0000, 


0. 


.9753, 


0. 


9512, 


0, 


.9277, 


0. 


.9048, 


0, 


,8825, 


0. 


8607, 


0. 


.8395, 


0. 


,8187, 


0. 




7985, 












































0. 


.7788, 


0. 


.7596, 


0. 


7408, 


0. 


.7225, 


0. 


.7047, 


0. 


.6873, 


0. 


6703, 


0. 


.6538, 


0, 


. 6376, 


0. 




6219, 












































0. 


.6065, 


0. 


.5916, 


0. 


5769, 


0. 


.5627, 


0. 


.5488, 


0. 


.5353, 


0. 


5220, 


0. 


.5092, 


0. 


,4966, 


0. 


\l 


4843, 












































0. 


.4724, 


0. 


.4607, 


0. 


4493, 


0. 


.4382, 


0, 


.4274, 


0. 


,4169, 


0. 


4066, 


0. 


.3965, 


0. 


,3867, 


0. 


\l 


3772, 












































0. 


.3679, 


0. 


.3588, 


0. 


3499, 


0. 


.3413, 


0, 


.3329, 


0. 


.3247, 


0. 


3166, 


0. 


.3088, 


0, 


,3012, 


0. 




2938, 












































0. 


.2865, 


0. 


.2794, 


0. 


2725, 


0. 


.2658, 


0. 


.2592, 


0. 


,2528, 


0. 


2466, 


0. 


.2405, 


0, 


,2346, 


0. 


)l 


2288, 












































0. 


.2231, 


0, 


,2176, 


0. 


2122, 


0. 


.2070, 


0, 


.2019, 


0. 


,1969, 


0. 


1920, 


0. 


.1873, 


0, 


,1827, 


0. 


\£ 


1782, 












































0. 


.1738, 


0. 


.1695, 


0, 


1653, 


0. 


.1612, 


0. 


.1572, 


0. 


.1534, 


0. 


1496, . 


0. 


.1459, 


0, 


,1423, 


0. 


kT 


1388, 












































0, 


.1353, 


0. 


.1320, 


0. 


1287, 


0. 


,1256, 


0. 


.1225, 


0. 


.1194, 


0. 


1165, 


0. 


.1136, 


0. 


,1108, 


0. 




1081, 












































0, 


.1054, 


0. 


.1028, 


0. 


1003, 


0. 


.0978, 


0, 


.0954, 


0. 


.0930, 


0. 


0907, 


0. 


.0885, 


0, 


,0863, 


0. 




0842, 












































0. 


.0821, 


0. 


.0801, 


0. 


0781, 


0. 


.0762, 


0, 


.0743, 


0. 


.0724, 


0. 


0707, 


0. 


.0689, 


0. 


.0672, 


0. 


^ 


0655, 












































0. 


.0639, 


0, 


.0623, 


0. 


0608, 


0. 


.0593, 


0. 


.0578, 


0. 


,0564, 


0. 


0550, 


0. 


.0537, 


0. 


.0523, 


0. 




0510, 












































0. 


.0498, 


0. 


.0486, 


0. 


0474, 


0. 


.0462, 


0, 


.0450, 


0, 


,0439, 


0. 


0429, 


0. 


.0418, 


0. 


.0408, 


0. 


)i 


0398, 












































0. 


.0388, 


0. 


.0378, 


0. 


0369, 


0. 


.0360, 


0, 


.0351, 


0. 


,0342, 


0. 


0334, 


0. 


.0325, 


0, 


.0317, 


0. 


1^ 



0310, 
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0. 


.0302, 


0. 


,0295, 


0. 


0287, 


0. 


.0280, 


0. 


,0273, 


0. 


.0266, 


0. 


.0260, 


0, 


.0253, 


0, 


.0247, 


0. 




0241, 












































0. 


.0235, 


0. 


.0229, 


0. 


0224, 


0. 


.0218, 


0. 


,0213, 


0. 


.0208, 


0. 


,0202, 


0, 


.0197, 


0, 


.0193, 


0. 




0188, 












































0- 


,0183, 


0, 


.0179, 


0. 


0174, 


0. 


.0170, 


0. 


.0166, 


0. 


.0162, 


0. 


,0158, 


0. 


.0154, 


0, 


.0150, 


0. 




0146, 












































0. 


,0143, 


0. 


.0139, 


0. 


0136, 


0. 


.0132, 


0. 


.0129, 


0. 


.0126, 


0. 


.0123, 


0. 


.0120, 


0, 


.0117, 


0. 




0114, 












































0. 


.0111, 


0. 


.0108, 


0. 


0106, 


0. 


.0103, 


0. 


,0101, 


0. 


.0098, 


0. 


.0096, 


0, 


.0093, 


0, 


.0091, 


0. 




0089, 












































0. 


.0087, 


0, 


,0084, 


0. 


0082, 


0, 


.0080, 


0. 


,0078, 


0. 


.0076, 


0. 


.0074, 


0. 


.0073, 


0, 


.0071, 


0. 




0069, 












































0. 


.0067, 


0. 


,0066, 


0. 


0064, 


0, 


.0063, 


0. 


,0061, 


0. 


.0059, 


0. 


.0058, 


0, 


.0057, 


0, 


.0055, 


0. 




0054 











































}; 

// compute the parameters for the color feature lookup table 

switch (pNbColor) { 

case kComplexityLow : 

lNbHueBin=NbColorBinBase ; 

lHueStepDW=255 . 0/ (double) INbHueBin; 

break; 

case kComplexityMedium: 

lNbHueBin=NbColorBinBase *2; 
lHueStepDW=255 . 0/ (double) INbHueBin; 

break; 



case kComplexityHigh : 

lNbHueBin=NbColorBinBase *4; 
lHueStepDW=255 . 0/ (double) INbHueBin; 

break; 

default: 

lNbHueBin=NbColorBinBase; 
lHueStepDW=255 . 0/ (double) INbHueBin; 

break; 



} 



lNbColFeature= lNbHueBin*2 +2; 

// resize the color feature array 

if (error = ResizeArraylD (pColSpectrumRef , sizeof (double ) , INbColFeature) ) 
{ 

goto END; 
} 

// Get the color array pointer 
if (error = GetArraylDPtr (pColSpectrumRef , &lColSpectrumAlDWPtr ) ) 
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{ 

goto END; 
} 

/*Initilization */ 
for (i=0; KlNbColFeature; { 
IColSpectrumAlDWPtr (i]=0; 

} 

// get the pointer to the input image structure 
ISrcIMPtr = pSrcIMRef; 
ISrcLineWidth = lSrcIMPtr->lineWidth; 

// if ROI input connected, calculate the color feature based on the ROI image 
if( pRoi->numContours>0) { 

if ( errors New Image ( &lMaskImagePtr , IMAGE_U8, 4, 4, 0) ) { 
goto END; 

} 

if (error = ROIToMask (IMasklmagePtr, NULL, pRoi, 1, &lpixelsOutsideMask) ) { 
goto END; 

}; 

IMskLineWidth = lMaskImagePtr->lineWidth; 
CalcMaskPart (pSrcIMRef , IMasklmagePtr , &lMaskPart) ; 

GetConstlmagePixel (pSrcIMRef , IMaskPart . startlmageX, IMaskPart . startlmageY, . & 
ISrcPixPtr) ; 

GetConstlmagePixel {IMasklmagePtr , IMaskPart . startMaskX, IMaskPart . startMaskY, wr 
(ConstPixelPtr*) &lMaskLinePtr) ; 

ISizeXLW = IMaskPart. sizeX; 
ISizeYLW = IMaskPart . sizeY; 
llmageAreaDW =0; /* init the Area */ 

// control negative sizes 

if (ISizeXLW < 0) 

ISizeXLW = 0; 

if (ISizeYLW < 0) 

ISizeYLW = 0; 

/* get image pixel base addresses */ 
// GetConstlmagePixel (pSrcIMRef , ICoord.xlLeft, ICoord.ylTop, &lSrcPixPtr) ; 
switch (lSrcIMPtr->imageType) { 
case IMAGE_RGB32: 

lSrcPixRGBPtr=lSrcPixPtr . PixRGB_Ptr ; 
break; 
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case IMAGE_HSL32: 

lSrcPixHSLPtr=lSrcPixPtr . PixHSL_Ptr ; 

breaks- 
default : 

error = ERR_BADIMAGETYPE; 

goto END; 



} 

// main switch 

switch (lSrcIMPtr->imageType) { 

case IMAGE_RGB32: 

for (ILineLW = 0; ILineLW < ISizeYLW; lLineLW++) { 

lBufSrcRGBPtr=lSrcPixRGBPtr; 
IMaskPixPtr = IMaskLinePtr ; 

for (IColLW = 0; IColLW < ISizeXLW; IC0ILW++) { 
if (*lMaskPixPtr++) { 
// initializatin 

lNuinColorEntry=0; 
// convert the RGB to HSL that corresponds to our board 



IRHHVal = lBufSrcRGBPtr->Byte_RGB.R 
IGSSVal = lBufSrcRGBPtr->Byte_RGB.G 
IBLVVal = lBufSrcRGBPtr->Byte_RGB.B 



CCONV_HSI {IRHHVal, IGSSVal, IBLVVal, Icoring, Ireplace, 1, s, h) 
Ihue = (Byte) (h +0.5 + offset) % 256; 
llgt = (Byte) (1+0.5) ; 
Isat = (Byte) (s + 0.5) ; 



if( lsat< BWHueThreshold ) { 

lNuinColorEntry=(llgt < BWLgtMiddle ) ? (lNbColFeature-2) : »^ 

(lNbColFeature-1) ; 

} 



else { 

// Threshold Black color using Nobelist's formular 

//lBlkThreshold= (long) ((128 - BlkLgtThreshold) *exp(-0.i^ 
025* (Isat - BWHueThreshold)) + BlkLgtThreshold ); 

ifdsat < 200 ) 

lBlkThreshold= (long) ((128 - BlkLgtThreshold) * ^ 
expLookup [Isat - BWHueThreshold] + BlkLgtThreshold ); 

else 

IBlkThreshold-BlkLgtThreshold; 
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if (llgt < IBlkThreshold ) { 

lNumColorEntry= ( lNbColFeature-2 ) ; 

} 

else { 

// offset the red color to the center of the first »^ 

bin 

Ihue = { Ihue + OffsetColor) % 255; 
// clustering non-black-white color feature 
lHueIndex= (double) Ihue / IHueStepDW; 

lNumColorEntry= {Isat < pSatThreshold) ? ( lHueIndex*2 ) : 

(lHueIndex*2+l) ; 

} 

} 

IColSpectruitiAlDWPtr [INumColorEntry] += 1.0; 
lImageAreaDW++; 
} // end of if (*lMaskPixPtr++ ) 

IBuf SrcRGBPtr++ ; 

} 

ISrcPixRGBPtr +=lSrcLineWidth; 
IMaskLinePtr += IMskLineWidth; 

} 

// calculate the discrete color spectrume 
for( i=0; i< INbColFeature; i++ ) { 

IColSpectrumAlDWPtr [i] /= llmageAreaDW; 

) 

break; 

case IMAGE_HSL32: 

for (ILineLW = 0; ILineLW < ISizeYLW; lLineLW++) { 

lBufSrcHSLPtr=lSrcPixHSLPtr; 
IMaskPixPtr - IMaskLinePtr; 

for {IColLW = 0; IColLW < ISizeXLW; IC0ILW++) { 
if {*lMaskPixPtr++) { 
// initialization 
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(lNbColFeature-1) ; 



lhue=lBufSrcHSLPtr->Byte_HSL.H; 
lsat=lBufSrcHSLPtr->Byte_HSL.S; 
llgt=lBufSrcHSLPtr->Byte_HSL.L; 

lNiimColorEntry=0; 

if{ lsat< BWHueThreshold ) { 

lNumColorEntry=(llgt < BWLgtMiddle ) ? (lNbColFeature-2) : 

} 

else { 

// Threshold Black color using Nobelist*s formular 

// lBlkThreshold= (long) ((128 - BlkLgtThreshold) *exp(-^ 
0.025*(lsat - BWHueThreshold)) + BlkLgtThreshold ); 

ifdsat < 200 ) 

lBlkThreshold= (long) ({128 - BlkLgtThreshold) * ^ 
expLookup [Isat - BWHueThreshold] + BlkLgtThreshold ); 

else 

lBlkThreshold=BlkLgtThreshold; 
if (llgt < IBlkThreshold ) { 

lNumColorEntry= ( lNbColFeature-2 ) ; 

} 

else { 

// offset the red color to the center of the \^ 



first bin 



(lHueIndex*2) : (lHueIndex*2+l ) ; 



} 



Ihue = ( Ihue + OffsetColor) % 255; 
// clustering non-black-white color feature 
lHueIndex= (double) Ihue / IHueStepDW; 
lNumColorEntry= (Isat < pSatThreshold) ? 



IColSpectrumAlDWPtr [INumColorEntry] += 1.0; 
lImageAreaDW++ ; 
} // end of if (*lMaskPixPtr++) 

lBufSrcHSLPtr++; 



lSrcPixHSLPtr+=lSrcLineWidth; 
IMaskLinePtr += IMskLineWidth; 
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} 

// calculate the discrete color spectrum 
for( i=0; i< INbColFeature; i++ ) { 

IColSpectruitiAlDWPtr [i] /= llmageAreaDW; 

} 

break; 

} // end of switch (lSrcIMPtr->imageType) 

} // end of ROI case 

else { // ROI not connected 

GetConstlmagePixel (pSrcIMRef , 0, 0, &lSrcPixPtr ) ; 

ISizeXLW = lSrcIMPtr->xRes; 
ISizeYLW = lSrcIMPtr->yRes; 

llmageAreaDW = (double) ISizeXLW * ISizeYLW; 

switch (lSrcIMPtr->imageType) { 

case IMAGE_RGB32: 

lSrcPixRGBPtr=lSrcPixPtr . PixRGB_Ptr ; 

break; 

case IMAGE_HSL32: 

lSrcPixHSLPtr==lSrcPixPtr . PixHSL_Ptr ; 

break; 

default: 

error = ERR_BADIMAGETYPE; 
goto END; 

) 

// main switch 

switch (lSrcIMPtr->imageType) { 
case IMAGE_RGB32: 

for (ILineLW = 0; ILineLW < ISizeYLW; lLineLW++) { 
IBuf SrcRGBPtr=lSrcPixRGBPt r ; 

for (IColLW = 0; IColLW < ISizeXLW; IC0ILW++) { 
// initializatin 
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lNumColorEntry=0; 
// convert the RGB to HSL that corresponds to IMAQ color board 

IRHHVal = lBufSrcRGBPtr->Byte_RGB.R 
IGSSVal = lBufSrcRGBPtr->Byte_RGB.G 
IBLVVai = lBufSrcRGBPtr->Byte_RGB.B 

CCONV_HSI (IRHHVal, IGSSVal , IBLVVal , Icoring, Ireplace, 1, s, h) 
Ihue = (Byte) (h +0.5 + offset) % 256; 
llgt = (Byte) (1+0.5) ; 
Isat = (Byte) (s + 0.5} ; 

if( lsat< BWHueThreshold ) { 

lNumColorEntry=(llgt < BWLgtMiddle ) ? (lNbColFeature-2) : \t 

(lNbColFeature-1) ; 

} 

else { 

// Threshold Black color using Nobelist's formular 

//lBl)cThreshold= (long) ((128 - BlkLgtThreshold) *exp(-0.wr 
025* (Isat - BWHueThreshold)) + BlkLgtThreshold ); 

if (Isat < 200 ) 

lBlkThreshold= (long) ((128 - BlkLgtThreshold) * ^ 
expLookup [Isat - BWHueThreshold] + BlkLgtThreshold ); 

else 

lBlkThreshold=BlkLgtThreshold; 
if (llgt < IBlkThreshold ) { 

lNumColorEntry= ( lNbColFeature-2 ) ; 

} 

else { 

// offset the red color to the center of the first ^ 



bin 



(lHueIndex*2+l) ; 



Ihue = ( Ihue + OffsetColor) % 255; 
// clustering non-black-white color feature 
lHueIndex= (double) Ihue / IHueStepDW; 

lNumColorEntry= (Isat < pSatThreshold) ? (lHueIndex*2 ) : y^ 



} 



} 



IColSpectrumAlDWPtr [INumColorEntry] += 1.0; 
lBufSrcRGBPtr++; 
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} 

ISrcPixRGBPtr +=lSrcLineWidth; 

} 

// calculate the discrete color spectruine 
for( i=0; i< INbColFeature; i++ ) { 

IColSpectrumAlDWPtr [i] /= llmageAreaDW; 

} 

break; 
case IMAGE_HSL32: 

for (ILineLW = 0; ILineLW < ISizeYLW; lLineLW++) { 
lBufSrcHSLPtr=lSrcPixHSLPtr; 

for (IColLW = 0; IColLW < ISizeXLW; IC0ILW++) { 

// initialization 

Ihue^lBuf SrcHSLPtr->By te_HSL . H 
lsat=lBuf SrcHSLPtr->Byte_HSL . S 
llgt=lBuf SrcHSLPtr->By te_HSL . L 

lNiiinColorEntry=0; 

if ( lsat< BWHueThreshold ) { 

lNumColorEntry=(llgt < BWLgtMiddle ) ? (lNbColFeature-2 ) : 

} 

else { 

// Threshold Black color using Nobelist's formular 

// lBlkThreshold= (long) ((128 - BlkLgtThreshold) *exp(-kr 
0.025*(lsat - BWHueThreshold)) + BlkLgtThreshold ); 

if (Isat < 200 ) 

lBlkThreshold= (long) ((128 - BlkLgtThreshold) * M 
expLookup[lsat - BWHueThreshold] + BlkLgtThreshold ); 

else 

lBlkThreshold=BlkLgtThreshold; 
if (llgt < IBlkThreshold ) { 

lNumColorEntry= (lNbColFeature-2 ) ; 

} 

else { 

// offset the red color to the center of the ^ 

first bin 



(lNbColFeature-1) 
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(lHueIndex*2) : ( lHueIndex*2+l ) 



Ihue = ( Ihue + OffsetColor) % 255; 
// clustering non-black-white color feature 
lHueIndex= (double) Ihue / IHueStepDW; 
lNuinColorEntry= (Isat < pSatThreshold) ? 



IColSpectrumAlDWPtr [INumColorEntry] += 1.0; 
lBufSrcHSLPtr++; 



} 

lSrcPixHSLPtr+=lSrcLineWidth; 

} 

// calculate the discrete color spectrum 
for( i=0; i< INbColFeature; i++ ) { 

IColSpectrumAlDWPtr [i] /= llmageAreaDW; 

} 

break; 

} // end of switch (lSrcIMPtr->imageType) 
} // end of if ( pRoi->numContours>0) 



END: 



return error; 

} 



